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Abstract. Two new equatorially antisymmetric solutions recently published by 
Ernst et al are studied. For both solutions the full set of metric functions is derived 
^ . in explicit analytic form and the behavior of the solutions on the symmetry axis is 

analyzed. It is shown in particular that two counter-rotating equal Kerr-Newman- 
NUT objects will be in equilibrium when the condition m? + = q'^ + b'^ is verified, 
W), whereas two counter-rotating equal masses endowed with arbitrary magnetic and 

electric dipole moments cannot reach equilibrium under any choice of the parameters, 
so that a massless strut between them will always be present. 
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In the recent papers [H [2] the notion of equatorial antisymmetry has been introduced for 
stationary axisymmetric electrovac spacetimes and two new exact solutions of that type 
have been constructed within the framework of SibgatuUin's method [3111]. Since Ernst et 
a/ presented their solutions only in terms of the Ernst complex potentials S and $ [5j, it 
would be of interest to have the complete metrics related to those solutions because they 
could facilitate the analysis of physical properties of the new equatorially antisymmetric 
spacetimes. Another aspect of the Ernst-Manko-Ruiz (EMR) solutions requiring the 
knowledge of the respective metric fields is the following: these solutions describe some 
binary systems of two counter-rotating masses in the presence of the electromagnetic 
field, so with the aid of the analytical expressions for the metric functions it would be 
possible to consider the equilibrium problems of two counter-rotating constituents in 
the EMR spacetimes. 

Bearing in mind the above two main objectives, in section 2 of the present paper 
the metrical fields of the EMR solutions will be derived using expansions of the 
determinantal formulae of paper [6] in the N = 2 case. Furthermore, in section 3 
the obtained analytical formulae fully defining the EMR spacetimes will be utilized for 
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the resolution of the equihbrium problem of two equal counter-rotating Kerr-Newman- 
NUT particles. Apart from the behavior of EMR solutions on the symmetry axis, the 
stationary limit surfaces, magnetic lines of force, ring singularities and some limits of 
these solutions will be also considered. 



2. The Ernst potentials and metric functions of EMR solutions 



The equatorially antisymmetric EMR solutions belong to the N = 2 subclass of the 
analytically extended multisoliton solution whose Ernst potentials S and $ in the general 
N = 2 case are given by the formulae [6] 
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where Pi are arbitrary complex parameters, a„ can take on arbitrary real values or occur 
in complex conjugate pairs, and a bar over a symbol means complex conjugation; hi{an) 
and f{an) are constant objects defined as 

hi{an) = ei + 2fif{an), /(a„) = f{z)\^=^^, (2) 

ci and /; being coefficients in the expressions of the potentials S and $ on the symmetry 
axis: 
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The first of the EMR solutions (henceforth referred to as solution I), representing 
two counter-rotating electrically and magnetically charged masses, is defined by the axis 
data 

, , z — k — m — iia + v) z + k — m + i(a — u) 



eiz 



k + m — i{a — z/) z + k + m + i{a + u) 
2{q + ib)z 



f(z) = "^'^ ^ (4) 

[z - k + m - i{a - u)][z + k + m + i{a + u)y 

where the parameters m, a, q, b are the mass, angular momentum per unit mass, electric 

charge and magnetic charge, respectively, u is the NUT parameter and k is half the 

coordinate distance between the masses. Using (jlj) one easily obtains the quantities 

2{m + iiy)[k — m + i(a — u)] 2{m + iu)[k + m + i(a + z/)] 
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I3i = k — m + i{a — u), (32 = —k — m — i(a + u), 

^ _ {q + ib)[k - m + i(a - z/)] j _ {l + i^) [k + m + i{a + v)] 

k + ia ' k + ia 

by decomposing e{z) and f{z) into simple fractions, and also the corresponding 
parameters namely, 

tti = a+, 02 = as = — 04 = — a+, 

a± = V5 ± 2c;, (5 = m2 + A;2^3z/2-a2_2(g2 + 62)^ 

= [(m^ + z/2 _ g2 _ ^2)(p ^ 2z/2 - - g2 _ _ _ fc^)2]i/2^ 

as roots of the algebraic equation 

e{z) + e{z)+2f{z)f{z) = Q. (7) 

Formulae ^ and ([HD fully determine the respective quantities and hi{an) 

that appear in the determinants ([1]). The potentials £ and $ calculated in [2] for the 
data (jl]) have the form 

^ = -. $ 



A + A + B' 

A= {m^ + -q^ - b^){[{m^ + P + + a^f - 4{mk + auf] 

x{R+ - i?_)(r+ - r_) - 5a+a.{R+ + /?_)(r+ + r_)} 

+2a+a_{[(m2 + z/^ - - b'^){m'^ - k'^ + - z/^) + 2(mz^ - ka)'^] 

x(_R+i?_ + r+r^) + 2id{mu — ka){R+R^ — r+r_)}, 

B = Ada^a-{m + iu){[m'^ + z/^ — — 6^ + i(mu — ka)] 

x(i?+ + i?_ - r+ - r^) - + + r+ + r_)}, 

C = 4(ia+Q;_(g + ib){[m'^ + z/^ — — 6^ + 2(mz/ — ka)] 

x{R+ + i?_ - r+ - r_) - d{R+ + i?_ + r+ + r_)}, 

i?± = v/p' + (z±a+)2, r± = + (z ± (8) 

The other EMR equatorially antisymmetric solution (henceforth solution II), 
describing a pair of counter-rotating masses endowed with electric and magnetic dipole 
moments, arises from the axis data 

z — k — m — i(a + u) z + k — m + i(a — u) 
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and the corresponding quantities 
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The Ernst potentials defined by the axis data iQ were found to have the form ^ 
A-B . C 



A = {6[{m^ + z/2)2 + (mz/ - /ta)^] - d\3m^ - P + + a^)} 

x{R+ - i?_)(r+ - r_) - a+a_[5(m2 + u^) - - c^] 

x(i?+ + R-){r+ + r_) + 2a+a_{[(m2 + z/^)^ + (mz/ - kaf - d^] 

X + r+r_) + 2id{mv - ka){R+R_ - r+r_)}, 

i? = Ada+a-{m + iz/){[m^ + z/^ + i(mz/ — ka)] 

x{R+ + R^ -r+ -r_) - + + r+ + r_)}, 

C = Ad{x + ic){[m^ + z^^ - i(mz/ - A;a)][a_(i?_ - - a+(r_ - r+)] 
+c/[a_(i?_ - + a+(r_ - r+)]}, 

i?± = v/p' + (^±a+)2, r± = + ± (12) 

The constants m, a, z/, in fl9|)- f|T2|) have the same meaning as in solution I, but c and 
X are the magnetic dipole and electric dipole parameters, respectively. 

The metric functions /, 7 and u which appear in the Papapetrou axisymmetric 
stationary line element 

ds^ = f-'[e^^dp^ + dz^) + p'^dif^] - f{dt - iodipf (13) 

can be calculated with the aid of the general formulae obtained in the paper [6]: 
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It is clear that in all practical applications of the above formulae for the metric 
functions, determinants ( |T5|) should be expanded and then evaluated for some particular 
axis data with the aid of a computer program for analytical calculations. In Appendix 
the reader can find the expansions of the determinants (ITSjl which have proved to be 
most efficient for symbolic computer processing. The results of the calculations of metric 
functions with the help of formulae (A.l), (A. 2) are given below. The functions /, 7 
and uj of both EMR solutions permit a unified representation, namely, 

AA-BB + CC o„ AA-BB + CC 



f 



.27 



{A + B){A + B)' 6Aala^„d^R+R_r+r. 
4iy + M^(^ + B)- G{A + B)- CI] 



(16) 



AA-BB + CC 

where the functions A, B, C, R±, r± and quantities a±, d are defined by formulae ([8] 
and (E]) in the case of solution I, and by formulae (fT2i) and (JTTl) in the case of solution II 
The functions G, H and /, entering the expression of u, have the following forn|j|: 

G = -zB + 2{m + w)A - 2da+a.{m^ + - - 6^) 

X [{a+ — a_)(-R_r_ — R+r+) + («+ + a_)(-R_r+ — -R+r_)] 
+4:da+a^[d — — z/^ + + 6^ + i{ka — mu)] 



H = zA 



X [a+{m + iz/)(i?_ — _R+) 
+4:da+a^ [d + + z/^ — 
X [a;_(m + iz/)(r_ — r+) - 
2(m + iz/)i? + a+a„(m — iu 



2{m^ + 
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b'^ — i(A;a — mu)] 



2{m^ + 
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x{(m2 + z/2 _ g2 _ [(^^ _ a_)^{R+r+ + R_r_) 

+ {a+ + a_)^(i?+r_ + i?_r+)] + 2d{m + iz/) 

X [(a+ — a_)(i?+r+ — R_rJ) + (a+ + a_)(i?+r_ — -R_r+)] 

+4(i?+i?_ + r+r_)[(m2 + _ ^2 _ ^2^)(^^2 _ ^2 ^ ^2 _ ^2^ 



-2(fca 



+ 8i(i(A;a — mv){Rj^R_ — r+r_)} 



+4(ia+Q;_(m + iz/){[m^ + z/^ — — 6^ — i(A;a — mu)] 
X The formulae of this paper have been obtained and checked using the Mathematica computer program 
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X — R_) — «-(r+ — r_)] — — RJ) + a;_(r+ — r_)]}, 

I = 2{q + ih){A + B)-[z + 2(m + w)\C + 2c/a+a_(m - iz/)(g + ih) 

X — a;_)(-R+r+ — R^r^) + (a+ + a„)(i?+r„ — -R_r+)] 
+4(ia+Q;_(g + i6){[m^ + z/^ — — 6^ — i(/ca — mu) — d\ 
x[2(m - w){R+ + + - - [m^ + _ ^2 _ ^2 

—i{ka — mv) + (i][2(m — iz^)(r_|_ + r_) + a-(r_|_ — r_)]}, (17) 
in the case of solution /, and 

G = -zB + 2(m + iz/)v4 + 2ci{(a+ - a_)[(m^ + z/^)a+a_ + + c'^]{R+r+ - R-T-) 
+ (a+ + a_)[(m^ + z/^)a+a_ - - c^](^+r_ - R^r+)} 
+4(ia;+(m + iz/){[(A; + ia)^ — (m + iz/)^][m^ + z/^ — c? + i{ka — mz/)] 
x(r_(_ — r_) — 2a_(m — iz/)[m^ + z/^ + d — i(A;a — mz/)](r_|. + r_)} 
+4:da^{m + iz/){[(A; + ia)^ — (m + iz/)^][m^ + z/^ + + i(/i;a — mz/)] 
x(i?_ - i?+) + 2a+(m - iz/)[m^ + z/^ - - i(A;a - mz/)](i?„ + i?+)}, 
= — 2(m + iz/)i? + (m — iz/)(a+ — «_){(«+ — a_)[(m^ + z/^)a_(.a_ + + c^] 
x(i?+i?_ + r+r„) + 2rfa;+a;_(m + iz/)(i?+i?_ — r+r.)} 
+ (m — iz/)(a+ + «_){(«+ + «_)[(m^ + z/^)a+a_ — X^ — c^] 
x(i?+r_ + R-r+) + 2(iQ;+a;_(m + iz^)(_R+r_ — i?_r+)} 
-Aa+a^{m - iz/){[(m^ + z/^)^ + {ka - mvf - dP]{R+R_ + r+r_) 
—2id{ka — mv){R^R_ — r^rS)} + 4da^a^{m + iz/) 
x{[m^ + z/^ - i(/i;a - mi^)][a_{r_ - r+) - a+(i?_ - i?+)] 
+c/[a_(r_ - r+) + a+(i?_ - i?+)]}, 

J = -{z + 2m + 2iz/)C + 2(m - iz/)(x + ic){[m^ + z/^ + i(A;a - mz/)][(a+a_ - 5) 

x(i?+r+ + R^r^) + (a+a_ + 5)(_R+r_ + R-r+) — 2«+a_ 

x{R+R. + r+r_)] + 2d2(i?+ - i?_)(r+ - r_) - 2da+a^{R+R- - r+r_)} 

+4(i(x + ic){a+a_[3(m^ + z/^) — d + i{ka — mv)]{r^ + r_) 

+2a+(m — iz/)[m^ + z/^ — (i + i(A;a — mz/)](r_|_ — r_) — 

x[3(m^ + z/2) + d + i{ka - mz/)](i?+ + i?_) - 2a_(m - iz/) 

x[m^ + z/^ + (i + i{ka - mu)]{R+ - R_) - 8da+a_{m + iz/)} (18) 

in the case of solution II. 

It is worthwhile mentioning that an arbitrary additive constant in the expression of 
u in ( |T6|) was chosen in such a way that the constant ujq in the definition of equatorially 
antisymmetric spacetimes [2] were equal to zero, i.e., u^p, z) = —uj{p, —z) automatically. 

3. The two— body equilibrium problem in EMR spacetimes 

The expressions of the metric functions obtained in the previous section can be used for 
the analysis of the equilibrium problem of two counter-rotating constituents of the Kerr- 
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Newman-NUT type. In order these constituents to be in equilibrium, it is necessary 
that the conditions P [IHl [HI [12] 

7 = 0, cu = (19) 
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Figure 1. Location of sources in the solutions I and II. 

are fulfilled on the part p = 0, \z\ < a_ of the symmetry axis that separates the particles 
(region II in Fig. 1). An important advantage of the equatorially antisymmetric systems 
is that the condition for u is fulfilled automatically in the region II, what can be checked 
directly with the aid of formulae f[T7|) and f[T8l) . so that the condition for the metric 
function 7 needs to be satisfied. The latter condition leads to the algebraic equation 
(m^ + z/2 - g2 _ 52^) p2 ^ ^2 ^ ^2 ^ ^2^)2 _ ^ _ ^^^_^] = 0, (20) 

in the case of solution I, and to the algebraic equations 

{5 - a+a_)[(m2 + + {ka - muf] + S[5 - A{m'^ + z/^) + (1 ± 2)a+a_] = 0, (21) 

in the case of solution II. It should be mentioned that for both solutions 7 = in the 
regions I and III, while uj = m the region I and uo = —Au in the region III due to 
the presence of two semi-infinite NUT singularities described in the paper [19] , the first 
singularity extending from a+ to +00 and the second from — a+ to —00. When z/ = 0, 
both conditions ([T9|l are verified identically in the regions I and III for solutions I and 
II because the latter solutions become asymptotically fiat in that limit. 
From (|20l) follows an important relation 

+ z/2 - g2 _ 52 ^ 0, (22) 

at which two counter-rotating Kerr-Newman-NUT constituents are in equilibrium 
independently of the distance between them. Formula f[22]) generalizes the balance 
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condition = which is verified by two Majumdar-Papapetrou equal charged masses 
[T3| [H] (condition = also determines two spinning charged masses of the Perjes- 
Israel-Wilson type [151 [161 [Hj). 

Furthermore, it can be shown that in the case of solution I there are no equilibrium 
states other than defined by the relation (l22l) because the equation 

(A;^ + + + u'^f - A{km + auf - = (23) 

has the roots that cause ci = in the denominator of the function 7 from ( JT71) and, 
therefore, should be discarded. Indeed, first putting the term to the right-hand 

side of equation ( l23l) and then taking square of the resulting equality, we arrive at the 
equation 

[{k - mf + (a - uf] [{k + mf + (a + uf] [{k^ + + a' + u^ f 

-4{km + az/)2 -{6 + 21)^ - Ab^ + 4b\6 + 26^)] = 0. (24) 
The first two factors of this equation become zero when 

k = m, a = V and k = —m, a = —u, (25) 

and in both cases the solution I reduces to the ordinary NUT solution electrically and 
magnetically charged (see, e.g., Ref. [H]), that is, the two-body problem degenerates 
to the case of a single body. Equating to zero the third factor and solving the resulting 
equation with respect to 6^, we get 

52 = i (k-^ + + 3z/2 - - 2g2 ± ^(p + ^2 + ^2 + ^^y _ ^f.^ + qj,)2^ _ (26) 

Substituting now (1261) in the expression of d defined by ([6]) it is straightforward to check 
that d = identically. 

In the case of solution II one can solve analytically the balance equations (l2Ti) in 
exactly the same way as equation (!20l) was solved for solution I. For each choice of the 
sign in fl2T|) one then comes to an algebraic equation of higher order than the initial 
one which, however, factorizes into three factors, one of which is d"^, and the other two 
permit, after equating them to zero, the solutions of the corresponding equations with 
respect to 6^ similar to formula ( l26l) . The direct substitution of the expressions for 6^ 
thus obtained into equations (pT!) shows that they do not satisfy the latter equations, 
thus being fictitious roots that must be discarded. The case d = must be discarded 
too on the same grounds as in the case of solution I. Therefore, in the equatorially 
antisymmetric systems of two magnetized masses there are no equilibrium states under 
any choice of parameters, and hence these masses are always supported by a massless 
strut [9] between them. 

It should be emphasized that solutions I and II are asymptotically fiat in the absence 
of the NUT parameter, in which case, by construction, they are regular on the upper 
and lower parts of the z-axis. In the presence of the NUT parameter, the conditions 
7 = 0, C(j 7^ are verified on those parts of the axis for both solutions, giving rise to 
two semi-infinite massive singularities of the NUT type [19]. We also mention that 
although it is tempting, following paper [20], to redefine the function 7, by adding a 
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specific constant, in such a way that the condition 7 = is verified on the intermediate 
part of the symmetry axis, this would only change the single supporting intermediate 
strut 7 7^ of finite extension to a pair of semi-infinite struts 7 7^ of Israel's type 
(massless in the absence of the NUT parameter), thus only worsening the situation and 
making the corresponding spacetimes frankly unphysical. 

4. The multipole moments, basic limits, stationary limit surfaces and ring 
singularities 

In the paper [2] the multipole structure of solutions I and II was not studied, so it would 
be of interest to clarify this characteristic of the EMR solutions in some detail. We 
have calculated the first four mass, angular momentum, electric and magnetic multipole 
moments (Mj, Jj, Qi and ifj, respectively) as these were defined by Simon ^Ij. During 
the calculations we have used the Hoenselaers-Perjes procedure [22] rectified by Sotiriou 
and Apostolatos [23], yielding the following expressions for the multipole moments: 

Mo = 2m, M2 = 2m(A;2 - + 3z/^ - a^) - Akau, Mi = M3 = 0; 

Jo = 2z/, J2 = 4fcma + 2z/(fc^ - 3m^ + z/^ - a^), = J3 = 0; 

go = 2g, Q2 = 2q{k'^-m'^ + u^~a^)-Ab{ka-miy), = Q-^ = 0; 

Ho = 2b, H2 = 2b{k'^ -m^ + u'^ - a^) +4:q{ka-mu), = = (27) 

(the case of solution T) and 

Mo = 2m, M2 = 2m(A;2 - m^ + Sz/^ - a^) - Akau, = M3 = 0; 

Jo = 2z/, J2 = 4fcma + 2z/(fc^ - 3m^ + z/^ - a^), Ji = J3 = 0; 

Qo = Q2 = 0, Qi = 2x, Q3 = 2xik'^ -m^ + u'^ -a^)-Ac{ka-mu); 

Hq = H2 = 0, Hi = 2c, H3 = 2c{P - m^ + - a^) + Ax{ka - mv) (28) 

(the case of solution IT). 

The above expressions support the physical meaning attributed to the parameters 
of EMR solutions in the paper |2J . From fl271) and fl28|l follows that the main difference 
between the two solutions lies in the structure of the electromagnetic moments: in 
solution I the odd moments Q2n+i and iJ2n+i are equal to zero, whereas in solution II 
are equal to zero the even moments Q2n and iJ2n- 

In the absence of the electromagnetic field, both EMR solutions reduce to the same 
special vacuum spacetime for two counter-rotating Kerr-NUT masses, and this limit 
belongs to the double-Kerr family of solutions of Kramer and Neugebauer [2^. It is 
interesting that by further setting k = m, a = u, one arrives at the single NUT solution 
[25] with the total mass 2m and the NUT parameter 2z/ which, as was demonstrated by 
Manko and Ruiz [19] , shares the property of being equatorially antisymmetric. Solution I 
reduces to the Breton-Manko electrovac solution [26] for two counter-rotating Kerr- 
Newman masses when the NUT parameter z/ and the magnetic charge b are equal to 
zero. The physically most interesting subclass of solution II is defined by vanishing 
parameters u and in which case this EMR solution represents two counter-rotating 
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Figure 2. Particular SLS of solution I demonstrating the growth of SLS with increasing 
NUT parameter v. In all three cases m = 4:, k = 3, a = l, q = b = 1/2, but v varies: 
(A) u = l,{B)u = 2, (C) u = 4. 



magnetized masses, giving the first example of solutions of this kind known in the 
literature. 
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Figure 3. Specific SLS of solution I: (A) The case of two hyperextreme constituents 
defined hy m = 2, k = a = 3, u = q = b = 1; (B) Two subextreme constituents with 
negative masses (m = —4, k = 3, a=l, i' = q = b = 1/2) and two ring singularities 
located at p ~ 1.2, z ~ ±5.935; (C) Two hyperextreme constituents with negative 

masses (m = —2, k = a = 3, i' = q = b= 1/2) and two ring singularities located at 
p ~ 3.246, z ~ ±3.801. The dots in (B) and (C) denote ring singularities. 



The stationary limit surfaces (SLS) which are defined by the equation / = 0, in 
the case of EMR spacetimes display several interesting properties worthy of mentioning 
here. In figures 2 and 3 we have plotted several particular SLS of solution I, while 
in figures 4 and 5 one can find particular SLS of solution IL Figures 2 and 4 clearly 
demonstrate a completely different evolution of SLS in solutions I and II as function 
of the NUT parameter: in solution I (see figure 2) the SLS grows with increasing u, 
while in solution II (figure 4) the increase of u causes the degeneration of SLS. In 
figures 3(A) and 5(A) the typical SLS of toroidal type show a clear similarity in the 
case of two hyperextreme constituents with positive masses for both solutions I and 
II. Nonetheless, in the case of constituents with negative masses the corresponding 
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ring singularities develop differently. Indeed, it follows from figure 3(B,C) that ring 
singularities accompanying the negative mass in solution I are located outside the SLS, 
though very close to it. At the same time, looking at figure 5(B,C), one can see that the 
ring singularities of solution II, either in the subextreme or in the hyperextreme cases, 
lie on the SLS, exactly as in the pure vacuum case. We remind that ring singularities 
arise as solutions of the equation A + B = 

(C) 

7.5 F ' ' ' ' ' ' 

5 
2.5 

N 

-2.5 
-5 

-7.5 [ 

0.5 1 1.5 2 2.5 3 
P 

Figure 4. Particular SLS of solution II demonstrating the degeneration of SLS with 
growing NUT parameter v. In all three cases to = 2, k — i, a — 1, c — x— 1/2, while 
V varies in the following way: (A) v = 1/2, (B) v ^1, (C) v = 3/2. 
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2 2.5 3 



7.5 
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2.5 


-2.5 
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(B) 



0.5 1 1.5 
P 



2 2.5 3 




0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 



P P P 

Figure 5. Specific SLS of solution II: (A) The case of two hyperextreme constituents 
defined hy m = 1, k = 3, a = 2, v = 1/4, c = x = 1/2; (B) Two subextreme 
constituents with negative masses (to = —4, k — 3, a = l, v — c — x~ 1/2) and two 
ring singularities located at p ~ 0.932, z ~ ±5.908; (C) Two hyperextreme constituents 
with negative masses (m = —2, k = a — 3, i> — c = x — 1/2) and two ring singularities 
located at p ~ 3.142, z ~ ±3.836. The dots in (B) and (C) denote ring singularities. 



5. Conclusions 

The new symmetry discovered and described by Ernst et al [1] permits a systematic 
study of a large class of counter-rotating masses within the framework of general 
relativity. Thanks to the paper pLj, and partly to the paper [19], we know for instance 
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that the well-known NUT solution belongs to the family of equatorially antisymmetric 
spacetimes; it is also clear now that the word 'antisymmetric' applied more than a 
decade ago by Breton and Manko to a system of two equal counter-rotating Kerr- 
Newman particles in [26j was quite appropriate. While the NUT solution represents a 
single body accompanied by two semi-infinite singularities, the new EMR 6-parameter 
solutions already describe the two-body systems endowed with singularities of the NUT 
type. In the present paper we have constructed all metrical fields for both EMR solutions 
and have solved analytically the associated equilibrium problems, obtaining the genuine 
equilibrium states only for a particular subclass of solution I defined by the relation 
+ z/^ = + 6^ (the two counter-rotating constituents become then hyperextreme). 
Some physical properties of the EMR spacetimes have been also studied. As a final 
remark we would like to observe that in view of the physical differences existing in the 
equatorially symmetric case between the systems with even and odd number of particles 
(see, e.g., [271 I2H1 [29]) it may be interesting to consider the three-body equatorially 
antisymmetric solutions and compare them with the EMR two-body spacetimes. 
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Appendix. Expansions of the determinants E±, F, G, H, I, Kq 

Since the determinants ([1]), (ITS!) contain the coordinates p and z only through functions 
r„, it is advantageous to expand these determinants over the two lines in which r„ appear, 
using the Laplace rule. The resulting expressions for the determinants employed in the 
analytical computer codes for obtaining formulae ffTT]) . f|T8l) have the form 

E± = A±T = AtB, 

A= 5^ {-iy^'rir,{ai-aj)RkRiRiRj 

l<t<j<4 

x[(afc - P2){ai - (3i)hi{ak)h2{ai) - {at - (3i){ai - (32)hi{ai)h2{ak)], 
{k<l; kjj^ij) 

4 

i=l 

x[(afc - P2){ai - (3i)hi{ak)h2{ai) - {at - /?i)(a« - /?2)/ii(a/)^2(«A:)] 
-Rk[{aj - (32){ai - (3i)hi{aj)h2{ai) - (a^ - (3i){ai - P2)hi{ai)h2{aj)] 
+Ri[{aj - (32){ak - /?i)/ii(«j)/i2(afc) - (ftj - (3i){oik - (^2)hi{ak)h2{aj)]}, 

4 

F = C = Y.i-'^y^'nRjRkRiMRjfiaj) 
1=1 
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x[{ak - p2){ai - Pi)hi{ak)h2{ai) - {ak - Pi){ai - p2)hi{ai)h2{ak)] 

-Rkf{ak)[{aj - P2){ai - Pi)hi{aj)h2{ai) 

-{aj - Pi){ai - P2)hi{ai)h2{aj)] 

+Rif{ai)[{aj - ^2){oik - Pi)hi{aj)h2{ak) 

-{aj - Pi){ak - P2)hi{ak)h2{aj)]}, 

(j <k< I; J, k,ly^i) 

l<i<j<4 

x[{ak - p2){ai - Pi)hi{ak)h2{ai) - {ak - Pi){ai - p2)hi{ai)h2{ak)] 

4 

+ J2i-'^y^'riRjRkRiRi{ajRj 
1=1 

x[{ak - P2){oii - Pi)hi{ak)h2{ai) - {ak - Pi){ai - P2)hi{ai)h2{ak)] 
-akRk[{aj - P2){ai - Pi)hi{aj)h2{ai) 
-{aj - P-i){ai - p2)hi{ai)h2{aj)] 
+aiRi[{aj - f32){ak - f3i)hi{aj)h2{ak) 
-{aj - (3i){ak - f32)hi{ak)h2{aj)]}, 

4 

H^zA-{Pi + f32)r + Y^{-lYr,aiRjRkRiRi 

1=1 

x{Rj[{ak - P2){ai - (3i)hi{ak)h2{ai) - {ak - I3i){ai - p2)hi{ai)h2{ak)] 
-Rk[{aj - P2){ai - (3i)hi{aj)h2{ai) - {aj - (3i){ai - p2)hi{ai)h2{aj)] 
+Ri[{aj - /?2)(afe - Pi)hi{aj)h2{ak) - {aj - Pi){ak - P2)hi{ak)h2{aj)]} 
- {-^y^'rirj{ai-aj)RkRiRiRj 

l<i<j<4 

x{Rk[ei{ai - Pi)h2{ai) - e2(a/ - /32)/ii(a/)] 
-Ri\ex{ak - (3i)h2{ak) - 62(0;*; - /32)hi{ak)]}, 
/=(/i + /2)(A-r) + (/?i + /32-^)F 

+ {-iy^'rirj{ai-aj)RkRiRiRj 

i<«<i<4 

x{e2[(ai - P2)hi{ai)Rkf{ak) - {ak - P2)hi{ak)Rif{ai))] 

-ei[{ai - Pi)h2{ai)Rkf{ak) - {ak - Pi)h2{ak)Rif{ai))]} 
4 

+ Yi-'^y^^riRjRkRiRi 

i=l 

x{[ei(Q;j - Pi)h2{aj) - 62(0^ - (32)hi{aj)][f{ak) - f{ai)]RkRi 
-[ei{ak - Pi)h2{ak) - e2{ak - (32)hi{ak)][f{aj) - f{ai)]RjRi 
+ [61(0;; - Pi)h2{ai) - 62(0;/ - p2)hi{ai)][f{aj) - f{ak)]RjRk} 
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4 

-iyr,aiRjRkRiRi{Rjf{aj) 

i=l 

- p2){o:i - Pi)hi{ak)h2{ai) - (a^ - - /?2)/ii(a0^2(afc)] 

-Rkf{ak)[{aj - /32){ai - Pi)hi{aj)h2{ai) 
-{aj - Pi){ai - P2)hi{ai)h2iaj)] 
+Rif{ai)[{aj - (32){ak - /9i)/ii(aj)/i2(«fc) 
-{aj - Pi){ak - P2)hi{ak)h2{aj)]} 

-R1R2R3R4 Yl 

l<i<j<4 

x[{ak - (32){ai - (3i)hi{ak)h2{ai) - {at - (3i){ai - P2)hi{ai)h2{ak)], 
Ko= Yl {-^Y^'{ai-aj)RkRlR^RA{o^k-B2){al-h)hl{ak)h2{a^ 

l<i<j<A 

-{oik - - I32)hi{ai)h2{ak)], (A.l) 

where 

i?fc = (afe-/3i)K-/52), Rk = {au - h){ak - (A.2) 

are the constant objects first introduced by Manko and Ruiz for the vacuum sohton 
solution [30j. Note that the expressions (A.l) are given without the common factor 

{Pl-P2)l{Yt=lRnRn). 
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